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INTRODUCTION 


It has been suggested that resonance absorption of low frequency 

magnetohydrodynamic waves could be responsible for heating the solar corona 
(lonson, 1982; Hollweg, 1984). In this paper an improved method for 
calculating the resonance absorption heating rate is discussed and the results 
are compared with observations in the solar corona. To accomplish this, the 
wave equation for a mildly dissipative, compressible plasma is derived from 
the linearized magnetohydrodynamic equations for a plasma with transverse 
Alfven speed gradients. For parameters representative of the solar corona, it 
is found that a two scale description of the wave motion is appropriate. The 
large scale motion, which can be approximated as nearly ideal, has a scale 
which is on the order of the width of the loop. The small scale wave, 
however, has a transverse scale much smaller than the width of the loop, 
0.3-250 km, and is highly dissipative. These two wave motions are coupled in 
a narrow resonance region in the loop where the global wave frequency equals 
the local Alfven wave frequency. Formally this coupling comes about from 

using the method of matched asymptotic expansions to match the inner and outer 
(small and large scale) solutions. The resultant heating rate can be 
calculated from either of these solutions. A formula derived using the outer 
(ideal) solution is presented, and shown to be consistent with observations of 
heating and line broadening in the solar corona. 

DERIVATION OF THE BASIC EQUATIONS 

The linearized momentum and induction equations can be written as 
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where a 2 = c 2 /47r a, /3 2 = rjl/p 0 and y 2 = iy|/p 0 . The functional forms for the 

viscosity coefficients, rj, and the electrical conductivity, a, are given in 
Braginskii (1965) . By neglecting derivatives of the dissipation coefficients 
themselves, it is implicitly assumed in the equations above that dissipation 
processes are important only in a narrow layer, and that over this layer terms 
which are proportional to derivatives of the dissipation coefficients can be 
neglected. This is easily verified a posteriori . It is further assumed that 
all variations in B 0 and /k are transverse to the field so that v A = v A (x) 
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only. The plasma is assumed to have fi « 1 so that the transverse gradients 
of B 0 can be neglected and the total pressure p can be approximated as p ~ 
Take the time derivative of the momentum equation and use the 
induction equation to eliminate all derivatives of B . Since dissipation is 
assumed to be small, one can neglect terms \hich are proportional to products 
of the dissipation coefficients. Assume that the loop is driven at z = 0 and 
that at z = L there is a perfectly reflecting boundary (see Hollweg, 1984 for 
more discussion on this point) . Fourier transform the equation using 


vj Jx,y,z,t) 
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This form automatically satisfies the boundary condition at z - L as long as 

kzn “ (2n + 1 )tt / (2L) with n - 0, 1, 2 Define the following 

dimensionless parameters V A 2 (x) = v A z (x)/v A0 , R = v A0 /wa 2 , Prll = j3 2 /a z , and 
Prl = 7 2 /a 2 , where v A0 is some typical value of the Alfven speed within the 
loop, R is the magnetic Reynolds number and Prl and Pr|| the magnetic Prandlt 
numbers associated with shear and compressional viscosity. With these 
definitions, all lengths measured in units of d the transverse scale of the 

loop, and /c 2 (x) = (v? d 2 /v A z (x) - l^ 2 ), the equation for the velocity can be 

written in the dimensionless form 
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The terms on the left hand side of this equation are the terms obtained from 
ideal MHD. Solutions of the ideal equation for various forms of the Alfven 
speed profile have been discussed by several authors in the context of solar 
coronal heating (Ionson, 1982; Rae and Roberts 1982) and in the context of 
solar wind acceleration in coronal holes by Davila (1985) . These 
investigations have demonstrated that this equation describes the propagation 
of MHD surface waves and guided wave modes in an inhomogeneous plasma. 

The terms on the right are all due to dissipation and are therefore non- ideal 
MHD terms. In the solar corona these terms are "small" in most locations 
since they are all multiplied by the inverse of the magnetic Reynolds number R 
which is on the order of 10 12 for typical solar parameters. The exception to 
this ordering is at the location x * x A where ac 2 (x a ) » 0. At this location 
the first term on the left can only be balanced by the dissipative terms on 
the right even though they are small. Sufficiently far from x A , in the outer 
region, the ideal MHD solution is a reasonable approximation to the actual 
solution. However, near ^ , in the inner region, the character of the 

solution changes (the equation changes from second to fourth order) and the 
assumption of ideal MHD is not valid. This is a classic example of a singular 
perturbation (Nayfeh, 1981) . In the following paragraphs singular 

perturbation theory and the method of matched asymptotic expansions will be 
used to obtain a solution which is valid both outside and inside the resonance 
layer. To do this one must first reduce the coupled equations for the vector 
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components of the velocity to a single equation for one velocity component, 
say v x . 

To illustrate the basic idea of the resonance absorption layer and to 
demonstrate the matched asymptotic solution method, let us assume that the 
dominant dissipation mechanism is shear viscosity. This assumption cannot be 
completely justified at this point so it must be regarded simply as an ansatz. 
An investigation which incorporates ohmic, compressive viscous and shear 
viscous dissipation self consistently in a numerical solution is currently 
underway. The results of this investigation will be published when they 
become available. Under these assumptions the wave equation can be written 

+ * Z Vi n = ie v Z v ±n ( 5 ) 


where e = (1 + Prl)/ (R V A 2 ). An equation for v xn correct to first order in e 
can be obtained for the apppropriate ordering 1 » ky 2 » k z 2 » eky 2 
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The outer solution is obtained by expanding the velocity, v xn , as a power 
series in the small parameter e. The lowest order term must then satisfy 
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Detailed solutions of this equation have been obtained before (Rae and 
Roberts, 1982; and Lee and Roberts, 1986). For our purposes it is only 
necessary to obtain the solution near the resonance layer, i.e. where x -► x A 
with x A defined by the resonance condition >c 2 (x A ) =0. In this region, the 
first term dominates and the solution can be approximated as 

v ^ ~ A In (x - x ; 
xn K a' 

where A is the wave amplitude determined by matching boundary conditions at 
the driver. 

The inner solution can be obtained by considering the scale stretching 
transformation given by f ** (x - )/a, where a is the small scale parameter 

to be determined in the problem. For the case presented here a 3 — e/X where 
^ = ~k z 2 d(ln p)/dx. Using the shear viscosity and resistivity coefficients 
given in Braginskii one can estimate a ~ 0.3 km. This is below the resolution 
of current instruments. Using the largest viscosity coefficient in Braginskii 
one can estimate km. If this second estimate is more nearly correct, 

these sheets could be observed with instruments with resolution on the order 
of 0.1 arc second. Using the transformation described above, the inner 

equation can be written 


447 


d 


( 9 ) 


r d ^ dv 

- - , - « - 
dc 1 dr J dc 


xn . . 2 2 

= - ik a Cv 

y xn 


Assume a power law expansion of v xn with a 2 as the expansion parameter, then 
the solution of the zeroth order equation is 
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In the limit as r 00 it can be shown that 


v (<;) ~ -iC In (x - x ) (11) 

xn 1 2 A 

This shows that by proper choice of the constants, namely C 2 = iA, the inner 
and outer solutions match as is physically required. 


Although in the steady state the details of the velocity profile within the 
resonance layer do depend on the dissipation mechanism (eqn 10) , the heating 
rate does not. The result presented below is obtained by considering the 
outer solution and integrating the Poynting flux over the entire surface of 
the resonace layer (Chen and Hasagawa, 1974) . An equivalent result could be 
obtained from the inner solution by integrating the volumetric heating rate of 
Braginskii (1965) over the volume of the resonance layer. The result is 


B , 1 d ( k ) ^ f . ^ ix 
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where x 1 and x 2 denote the positions of the two surfaces of the resonance 
layer and A s is the surface area of the resonance layer. If we consider the 
limit of this equation as (x 1 -X 2 ) -* 0 and use causality to choose the proper 
analytic continuation (McPherson and Pridmore- Brown, 1966; Kapraff and 
Tataronis , 1977; Mok and Einaudi , 1985; Einaudi and Mbk, 1985; Bertin et al, 
1986) the resulting heating rate is 

2 ^2 

H = A — — x \a\ Z (13) 

s 8 mod k 2 

7 

To determine whether this heating rate is consistent with observational 
constraints, let us simply equate it to the observed radiation rate in soft 
x-rays to obtain an expression for the RMS velocity amplitude required to 
explain the observed emission. Using typical values of the parameters, B 0 ~100 
G, d=5xl0 8 cm, A(T)=1(T 22 (Rosner et al . 1978), P=300 sec, k y =2k z and Jn e 2 dl = 
10 2 8 * 2 9 (Webb et al. 1986) one obtains an estimate of V rms = 2-6 km/sec using 
the following expression. 
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Thi s is comparable to the observed value of 10-20 km/sec regularly seen from 
observations of non- thermal line broadening in the corona. 

CONCLUSIONS 

The primary conclusion to be drawn from these calculations is that to the 
level of the approximation adopted here, the observations of the heating rate 
and non- thermal line broadening in the solar corona are consistent with 
heating by the resonance absorption mechanism. This basic agreement is 
gratifying but several problems remain. The plane symmetry assumed here is 
highly idealized. It has been shown in the plasma physics literature that for 
the tokamak problem introducing cylindrical symmetry has presented no new 
physics. Nevertheless, when considering the heating rate to an accuracy of 
say factors of 2-5 the geometry factors must be properly accounted for. 
Second, observations of the turbulent power spectrum at the base of the corona 
are badly needed as input for the theory. These observations should be 
carried out in ions which are present at or above the transition region 
temperatures. It seems that EUV observations would be the most appropriate. 
Finally, other sources of dissipation must be considered. For although the 
heating rate is independent of the dissipation mechanism for any reasonable 
value of the coefficients, the amplitude of the velocity inside the resonance 
layer, and the width of the layer, both depend on the magnitude of the 

dissipation coefficient. High resolution instruments such as POF may be able 
to observe velocities within the narrow resonace regions in the reasonably 

near future. Therefore it is worthwhile to consider theoretically the 

observational consequences of various dissipation mechanisms now. In the 
work presented here, shear viscosity was assumed to be the dominant 

dissipation mechanism. This is not necessarily the case for solar 
conditions . 
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